Journal of Research on the Lepidoptera 12(3):135-144, 1973 


1160 W. Orange Grove Ave., Arcadia, California 91006, U.S.A. 
© Copyright 1973 


THE STUDY OF 
POPULATIONS OF LEPIDOPTERA 
BY CAPTURE-RECAPTURE METHODS 


P. M. SHEPPARD and J. A. BISHOP 
Department of Genetics and Department of Zoology, 
University of Liverpool, Liverpool L69 3BX, England. 
INTRODUCTION 

REcENTLY Kirpurr (1973) presented mark-release-recapture 
data on a population of the butterfly Euptychia hermes from 
which he determined daily population size and a maximum 
life expectancy. Fosdick (1973) has used identical methods 
in his study of the butterfly Anartia amalthea. Unfortunately 
those employed for analysing the population size are inefficient 
when such extensive multiple recapture data are available, and 
in addition they violate the assumptions implicit in the under- 
lying mathematical model. Furthermore, the method used for 
estimating the maximum life expectancy is spurious, since it 
ignores variation in the number of adults released and the 
ease with which butterflies were captured from day to day. 

Since many people are now using mark-release-recapture 
methods for studying populations we think it appropriate to 
mention some of the best available modern methods in order 
that people can analyse their data efficiently. For the purpose 
of this paper it would have been convenient to analyse the 
data of Kilduff and Fosdick but, unfortunately, they are given 
in such a way that they cannot be analysed by modern methods. 
Thus, Kilduff tells us that the numbers in the body of his trellis 
diagram refer to individuals recaptured. However, where a 
butterfly has been recaptured more than once it is not clear 
whether the entry in the diamond is appropriate for the first 
time the insect was captured or the last time, or if recaptured 
more than twice for some intermediate time. Fosdick declares 
that in his diagram the numbers also refer to individua!s. Never- 
theless from Table 1 in his paper it can be seen that the total 
(189) is greater than the total number of individuals (92). 
Consequently it seems appropriate not only to suggest efficient 
methods of analysis but also how the data should be presented 
in order that the reader may analyse it by modern methods 
should the author not care to do so. 
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CAPTURE-RECAPTURE METHODS 


Many populations of mobile animals can be investigated 
using capture-recapture methods. The intuitively easy to under- 
stand Lincoln Index (or Petersen estimator) is the simplest and 
is the basis of more complex mathematical models 


le = "YE Ie 
Rp 
where 
P is the population estimate on days 1 and/or 2 (see 


below) 


M, is the number marked and released into the pop- 
ulation on day 1 


N, is the number captured on day 2 


P 


Rə is the number of recaptures on day 2 
When R; is large (i.e. over about 20) the estimate of P is fairly 
unbiased but if there are fewer recaptures Bailey (1951, 1952) 
has shown that the estimate tends to be too large and is less 
biased if 1 is added to N, and R.. 


PE MN NAAN 
(R:+1) 
This correction can also be used in the modified Lincoln Index 
incorporated in multiple recapture models. 


The Lincoln Index makes certain assumptions: 
1) Marking does not affect animals and marks are not lost. 
2) The marked animals mix randomly in the population. 


3) The population is sampled at random. An animal's age, 
sex or mark does not influence its chance of being cap- 
tured (see later). 


4) Either there is no birth or immigration (dilution) or 
there is no death or emigration (loss). If there is no 
dilution but loss the estimate is valid for the population 
on the day of release. On the other hand if there is 
dilution but no death then the estimate is valid for the 
day of recapture. 


0) Sampling must be performed at discrete time intervals 
and the time spent sampling must be negligible (in real 
terms this means that it must be small) in relation to 
the inter-sampling period. 
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Most of the multiple capture-recapture models derived from 
the Lincoln Index try to remove assumption 4 to make allowance 
for dilution and loss and they more obviously attempt to repre- 
sent a real animal population. The most popular approaches are 
listed in Table 1. In general, as one proceeds down this table, 
one is dealing with methods that need increasing amounts of 
data (in the form of recaptures). However the models are in- 
creasingly realistic and give a more complete series of estimates 
of parameters of the population. In a capture-recapture study, 
the type of analysis will depend on the amount of data that can 
be obtained. Thus where comparatively few recaptures are 
made Fisher and Ford's method may be the most appropriate, 
but where more data are available Jolly's, or even Manly & Parr's 
method may be used. 


It is important to realise the limitations imposed by assump- 
tion 3 (above). If all individuals are not equally catchable vast 
errors can occur (see Bishop and Bradley 1973). Any data for 
analysis must be taken from a population that is as homogeneous 
as possible. This may demand separate analyses of sub-popula- 
tions (e.g. The sexes may have to be treated independently.) 
Any one sub-group may show variability (e.g. polymorphism) 
that will influence the probability of capture or recapture, so 
further divisions may be needed. Zoological insight used to 
assemble data from organisms into homogeneous groups will im- 
prove estimates of all parameters obtained from the population 
models listed in Table 1. If the nature and limitations of a 
method are understood then it is frequently convenient to use 
a computer to analyse data. Sources of suitable computer pro- 
grams are tabulated. 


EXAMPLES AND DISCUSSION 


Dowdeswell, Fisher and Ford (1940) provide data from a 
population of the blue butterfly Polyommatus icarus on Tean, 
Isles of Scilly, U.K. These have been presented in such a way 
that they can be re-analysed to show the sort of estimates that 
can be obtained from Fisher and Ford’s (1947) model (Table 2). 


The units in the Fisher and Ford data trellis (Table 2) are 
marks not animals (it is a Type A trellis in the sense of Leslie 
and Chitty (1951)). A butterfly caught on 26th August and 
recaptured on 29th and 30th August contribues one “mark” to 
29th August and two “marks” to 30th August. This model also 
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Table 2 


Fisher and Ford type data trellis and estimates of parameters for a population of 
the blue butterfly Polyommatus icarus on Tean, Isles of Scilly, U.K., 1938. (Data 
from Dowdeswell, Fisher and Ford 1940) Bailey's correction has been used in the 


calculation of population estimates. 





Aug. 
Captures Releases 26 27 3 6 7 
26 Aug. 1938 40 ho 
27 kz RO 5 
29 13 42 O 3 
30 52 50 3 8 5 
31 56 51 6 42 6 15 
2 Sept. 52 52 hà — 10 3 16 . 4} 
3 50 50 4 5 ae 44 5 qb 
6 15 15 1 1 1 5 4 5 5 
7 20 20 4 1 2 3 2 7 8 6 
7 e 0 O o 2 2 4 4 


N.B. The population was not sampled on 28th Auge, 1st, 4th, 5th Sept. 


TOTAL DAYS SURVIVED BY MARKS = 786 
DAILY SURVIVAL RATE = — 0.82851 


27 Aug 25 13.6 
28 61 35 
29 50 28.6 33 
30 52 32.1 29 
34 84 70.2 25 
1 Sept 112 21 
e 95 90.6 48 
te 120 80.4 20 
| ^ 4h4 21 
5 117 19 
6 97 112.5 16 
| 7 92 147.6 15 
| 8 





TOTAL GAINS = 116 
TOTAL LOSSES =262 


* Third order harmonic running means. + Calculated using smoothed population 
estimates. 
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Table 5 


Basic capture-recapture data collected from a population of the burnet moth Zygaena 
filipendulae at Dale, Pembrokeshire, U.K. The complete history of each recaptured 
individual is shown (capture or recapture indicated by 1) (from Manly & Parr 1968). 


JULY 1968 JULY 1968 
19 20 21 22 24 19 20 21 22 24 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 il 1 
1 1 1 1 1 1 
1 2 1 1 1 1 
1 1 1 1 1 
1 1 1 1 1 
1 1 1 1 1 
1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 
1 1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 NEW E 1 1 
1 1 1 1 1 1 
1 1 1 1 1 1 
1 1 1 1 4 1 
1 1 1 1 
1 1 1 1 1 
1 1 
Number captured once only 
(on day indicated) and not 2 qo 5 9 1 


included above. 
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Table 4 


Type À and Type B data trellises derived from data in Table 5. 


Type A data trellis (for Fisher and Ford's model. Units under 'Date of Marking! 


refer to marks not animals. 


Date of marking 


Date Captured Released 19" 20 FAI 22 
19 57 57 
20 52 52 22 
21 52 52 19 20 
22 51 31 9 6 142 
24 54 54 w B æ 15 


Type B data trellis (for Jolly's model). Units record date of last capture of 


an animal recaptured. 


Date of marking 


Date Captured Released Ho A 21 ee 
19 en 27 
20 52 52 25 
21 52 52 8 20 
22 31 31 2 29 We 


24 54 54 1 6 WEM B 
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Table 2 


Estimates of population size and daily survival rate obtained from data of Table 5 
by Fisher and Ford's, Jolly's and Manly and Parr's models. (Bailey's correction 


has not been used) 


pate BHEEUES 6 mE JOLLY MANLY AND PARR 

July 1968 Population Size Pop.Size ae Bep Se peniya] 

19 E = 0.78 | - 0.76 

20 91 95 0.74 90 0.71 

21 102 98 0.76 97 0.69 

22 116 106 z 93 " 

23 2 - - 3 i 

eh 75 = » 


Fisher & Ford survival 
= 0.801 


The expectation of life is approximately 43 days. This can be estimated by 


1 ; s : 21 
deSürvivAm cate if all deaths occur just before sampling or by log. DOPVAVAMERESS 


if deaths are at random throughout the intersampling period. 


requires a quantity called the total days survived by all marks. 
On the 29th the mark of the 26th has survived 3 days and on 
the 30th 4 days. In addition the new mark of the 29th has sur- 
vived one day. Therefore this animal contributes 3 + 4 + 1 days 
to the total. The units in the Jolly data trellis refer to individual 
animals and the date they were last captured (i.e. it is a Type 
B trellis). Examples of Types A and B are given in Table 4. 
The Fisher and Ford method uses the information in the data 
trellis to calculate a 'survival rate per unit time for marked 
animals. This can then be used to allow for loss of marks (by 
death or emigration) in a modified Lincoln Index calculation. 
Estimates of population size then can be obtained and gains 
(birth and immigration) calculated. The complete analysis is 
given in Table 2. A worked example of Dowdeswell, Fisher and 
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Ford's data by Mr. R. J. White, Department of Zoology, Uni- 
versity of Liverpool, is available from J. A. Bishop. The calcu- 
lations for Jolly’s method are given by Jolly (1965) and South- 
wood (1966). 


Manly and Parr (1968) present a set of results from a cap- 
ture-recapture investigation of the burnet moth Zygaena filipen- 
dulae. They also provide a worked example of their model 
based on this data. Their data are tabulated in a manner that 
allows re-analysis by Fisher and Ford's and by Jolly's methods. 
The original data and the Type A and Type B trellises appear 
in Tables 3 & 4. The estimates of the population parameters 
provided by all three models are also given (Table 5). 


We believe that these examples emphasize the advantages 
of modern capture-recapture methods. They also indicate the 
desirability of presenting, where possible, the original data in 
a simplified but coherent form ( Table 3). 
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